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^/^ ■ Abstract 

OO ■ In a similar fashion to estimates shown for Harmonic, Wachspress, and Sibson coordinates 

in [Gillette et al., AiCM, to appear], we prove interpolation error estimates for the mean value 
coordinates on convex polygons suitable for standard finite element analysis. Our analysis is 
, based on providing a uniform bound on the gradient of the mean value functions for all convex 

polygons of diameter one satisfying certain simple geometric restrictions. This work makes 
rigorous an observed practical advantage of the mean value coordinates: unlike Wachspress 
' coordinates, the gradients of the mean value coordinates do not become large as interior angles 

, of the polygon approach tt. 
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1 Introduction 

Barycentric coordinates are a fundamental tool for a wide variety of applications employing trian- 
00 ■ gular meshes. In addition to providing a basis for the linear finite element, barycentric coordinates 

I also underlie the definition of higher-order basis functions, the Bezier triangle in computer aided- 

• design, and many interpolation and shading techniques in computer graphics. The versatility of 

this construction has led to research attempting to extend the key properties of barycentric coordi- 
nates to more general shapes; the resulting functions are called generalized barycentric coordinates 
(GBCs). Barycentric coordinates are unique over triangles [37], but many different GBCs exist 
for polygons with four or more sides. While GBCs have been constructed for non-convex poly- 

■ gons (UllHlES] and smooth shapes [71 [HI [Ml El] , the most complete theory and largest number of 

■ GBCs exist for convex polygons. 
Interpolation properties of barycentric coordinates over triangles have been fully character- 
ized |22| 117]. Interpolation using GBCs, however, has a more complex dependence on polygonal 
geometry. The earliest GBC construction, now called the Wachspress coordinates [36], exhibits the 
subtleties of geometrical dependence: if the polygon contains interior angles near tt, gradients of 
the coordinates become very large. The more modern mean value coordinates [11] seem to avoid 
this problem. Floater et al. exhibit a series of numerical experiments showing good behavior of 
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Figure 1: Comparison of Wachspress and mean value coordinates over two pentagons with vertices 
(—1,1), (—1,-1), (1,-1), (1,1), and (0, x) where the value x is indicated in the figure. The 
coordinate for the final vertex (0, x) is plotted. For x = 1.5, the polygon contains no large interior 
angles and the gradient of both coordinates is well-behaved. As x approaches 1, the interior angle 
at (0, x) approaches vr and the Wachspress coordinate becomes very steep while the mean value 
coordinate has a bounded gradient. 

the gradients of mean value coordinates on polygons with interior angles close to vr [12]. The dif- 
ference in behavior can be observed on a very simple polygon as shown in Figure [H Combining 
well-behaved gradients with a simple and explicit formula, the mean value coordinates have be- 
come quite popular in the computer graphics community [IHl [23l [U [301 Ell [28] . Additionally, they 
have been implemented in finite element systems where they produce optimal convergence rates in 
numerical experiments [331 [Ml [39] . 

Our aim in this work is to mathematically justify the experimentally observed properties of mean 
value coordinates by proving a bound on their gradients in terms of geometrical properties of the 
polygonal domain. The gradient bound allows us to prove an optimal order error estimate for finite 
element methods employing mean value coordinates over planar polygonal meshes. Our approach 
follows that of our previous work [16] . where we carried out a similar program for other types of 
GBCs previously proposed for use in the finite element context: Wachspress [36], Sibson [HH [33l 
[32l I27j . and Harmonic coordinates [201 126j . Note that gradients of a ID rational interpolant with 
certain similarity to the mean value coordinates have been shown in [141 [T]. but gradients of the 
mean value coordinates themselves have not been analyzed previously. 

Our error estimate is contingent upon the mesh satisfying two geometric quality bounds: a 
maximum bound on element aspect ratio and a minimum bound on the length of any element edge. 
These are the same hypotheses assumed for our prior analysis of Sibson coordinates, placing the 
two coordinate types on par with regard to convergence in Sobolev norms. For scientific computing 
purposes, however, the mean value coordinates have several advantages. While Sibson coordinates 
are only continuous on the interior of an element |311 llOj. the mean value coordinates are C°°, 
reducing the complexity of numerical quadrature schemes required for their use. Sibson coordinates 
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Figure 2: Notation used in the paper. Vertices are always denoted in boldface. 



also require the construction of the Voronoi diagram while mean value coordinates are defined by 
an explicit rational function. This straightforward definition also allows mean value coordinates to 
be computed for non-convex polygons [H]; we comment on the applicability of our analysis in the 
non-convex setting in the conclusion. 

The remainder of the paper is organized as follows. In Section [21 we fix notation and review 
relevant background on polygonal geometry, mean value coordinates, and interpolation theory in 
Sobolev spaces. In Section [3l we establish a number of initial estimates on various quantities 
appearing in the definition of the mean value coordinates. Our main result is Theorem 14.31 in 
Section H] which gives a constant bound on the gradients of the mean value coordinates given two 
specific geometric hypotheses. As established in Lemma 12.71 this bound suffices to ensure the 
desired optimal convergence estimate, even when interior angles are close to vr. We give a simple 
numerical example and discuss applications of our analysis in Section [5l 

2 Background 

2.1 Polygonal Geometry 

Mean value coordinates will be analyzed in the same setting as |16j . We briefiy outline the primary 
notation and definitions. Let Q be a generic convex polygon in R^. The n vertices of Vt are denoted 
by vi, V2, . . ., v„ and let the interior angle at Vj be j3i] see Figure El The largest distance between 
two points in $7 (the diameter of is denoted diam(r2) and the radius of the largest inscribed 
circle is denoted p(f^). The aspect ratio (or chunkiness parameter) 7 is the ratio of the diameter 
to the radius of the largest inscribed circle, i.e., 7 := d\a,m.{Q) / p{Q) . 

Interpolation error estimates involve constraints on polygon geometry. For triangles, the most 
common restrictions bound the triangle aspect ratio or exclude triangles with angles smaller /larger 
than a given threshold. Generalizing this idea to convex polygons leads to a richer collection of 
potential geometric constrains, as many are no longer equivalent. For example, a bound on polygon 
aspect ratio does not imply an upper bound on interior angles. The two geometric constraints that 
we will require for establishing error estimates are listed below. 

Gl. Bounded aspect ratio: There exists 7* G M such that 7 < 7*. 

G2. Minimum edge length: There exists (i^, G M such that |vj — Vj| > > for all i ^ j. 
A third constraint restricting the maximum interior angle was used in |16j . 
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G3. Maximum interior angle: There exists /3* S M such that (3i < P* < tt for all i. 



While Gl3]was necessary in the analysis of Wachspress coordinates, we emphasize that Gl3]is not 
used in our present analysis of mean value coordinates. In fact, insensitivity to large interior angles 
is one of the primary motivations for using mean value coordinates |15] . By establishing an error 
estimate without assuming Q3] gives a stronger theoretical justification for this original motivation 
. In |16] we showed that under Gl and G2, two other closely related properties also hold. 

G4. Minimum interior angle: There exists G M such that /3j > /3* > for all i. 

G5. Maximum vertex count: There exists n* G M such that n < n*. 

Proposition 2.1 (Proposition 4 in [IB]). Under Gl and G2, 04 and 05 hold as well. 

Hence, when assuming only Gl and G2 for our analysis, we may also use G4 and G5 if needed. 
Assuming GUJand G|2l a sufficiently small ball cannot intersect two non-adjacent segments as stated 
precisely in the following proposition. 

Proposition 2.2 (Proposition 9 in [16j). There exists /i* > such that for all unit diameter, 
convex polygons satisfying CJH an-d (JH and for all x £ Q, i?(x, /i^,) does not intersect any three 
edges or any two non-adjacent edges ofVL. 

Remark 2.3. Restricting to only diameter one polygons in Proposition \2.S\ is sufficient for our 
analyses due to the (forthcoming) invariance property iJ3 

2.2 Generalized Barycentric Coordinates 

Barycentric coordinates on general polygons are any set of functions satisfying certain key properties 
of the regular barycentric functions for triangles. 

Definition 2.4. Functions Aj : — )■ M, i = 1, . . . , n are barycentric coordinates on Q. if they 
satisfy two properties. 

Bl. Non-negative: Aj > on 0. 



Most commonly used barycentric coordinates, including the mean value coordinates, are invari- 
ant under rigid transformation and simple scaling which we state precisely. Let T : — )• be 
a composition of rotation, translation, and uniform scaling transformations and let {A^} denote a 
set of barycentric coordinates on TVL. 

B3. Invariance: Ai(x) = Af (T(x)). 

Remark 2.5. The invariance property can he easily passed through Soholev norms and semi-norms, 
allowing attention to he restricted to domains $7 with diameter one without loss of generality. The 
essential case in our analysis is the H^-norm (defined more generally in Section \2.4\ ), l^lHi(rj) ~ 



n 



B2. Linear Completeness: For any linear function L : — )■ M, L 





where Vu = {du/dx,du/dy)'^ is the vector of first partial derivatives of u, and 
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for simplicity T is a uniform transformation, T(x) := /ix. Throughout our analysis, the Euclidean 
norm of vectors will be denoted with single bars \-\ without any subscript. Applying the chain rule 
and change of variables in the integral gives the equality: 



[ |VAf(x)|'dx= / 

h"-^ [ \VX^{y)\'dy=h''~^\X. 
Jn 



^V(Af(/.x))" 



dx 



The scaling factor h'^ resulting from the Jacobian when changing variables in the integral is the 
same for any Sobolev seminorm, while the factor of from the chain rule depends on the order of 
differentiation in the norm (1, in this case) and the semi- norm used {p = 2, in this case). When 
developing interpolation error estimates, which are ratios of Sobolev norms, the former term (i.e., 
the chain of variables portion) cancels out and latter term (i.e., the chain rule portion) determines 
the convergence rate. 

Several other familiar properties immediately result from the definition of generalized barycen- 
tric coordinates (Bl and B2). 

n 

B4. Partition of unity: Aj = 1. 

i=l 
n 

B5. Linear precision: ^^VjAj(x) = x. 

i=l 

B6. Interpolation: Aj(vj) = 5ij. 

Proposition 2.6. Suppose Bl and B2 hold. Then B4, B5, and B6 hold as well. 

Proof. B4 and B5 are merely special cases of B2, viz. L(x) = 1 and L(x) = x, respectively. 
Substituting x = Vj into B5 gives VjAi(vj) = Vj. Since is convex, this equality can only hold 
when Aj(vj) = 8ij. □ 

Having outlined the generic properties of generalized barycentric coordinates, we can now turn 
to the specific construction in question. 

2.3 Mean Value Coordinates 

The mean value coordinates were introduced by Floater [llj (see also |12| and the 3D exten- 
sion [l5])- For a point x in the interior of fi, define angles ai(x) := ZvjXVj+i and distances 
rj(x) := |x — Vj|; see Figure [2j Then for vertex Vj, a weight function Wi{'x.) is given by 



'"(^)+""'(^) ^ t,(x)+i._.(x) 
|vi-x| ■ ri(x) 



where ti(x) := tan(aj(x)/2) is used to simplify the notation. The mean value coordinates are given 
by the relative ratio of weight functions of the different vertices: 
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As in |16| . the primary task in developing interpolation estimates for a particular coordinate is 
bounding the gradient of the coordinate functions. The primary challenge with mean value coor- 
dinates stems from the fact that the weight functions wt are unbounded over the domain; when 
rj(x) approaches zero near vertex Vj or ai(x) approaches tt near boundary segment VjVj+i, Wi can 
be arbitrarily large. As we show in Theorem 14.31 however, this behavior is always balanced by the 
summation of weight functions in the denominator of Aj, resulting in a bounded gradient. 

2.4 Interpolation in Sobolev Spaces 

We set out notation for multivariate calculus: for multi-index a = (01,02) and point x = 
define x" := x^^y^^, a! := aia2, \a\ '■= ai +0^2, and D°'u := d^°'^u/dx°'^dy°'^ . In this notation, the 
gradient, i.e. the vector of first partial derivatives, can be expressed by 



Z)(o,i)^ 



The Sobolev semi-norms and norms over an open set are defined by 



I^l//'"(r2) / ^ |Z)"'u(x)pdx and ^ I'"l/f™(r2) • 

-'^ \a\=m 0<k<m 

The H^-norm is the L^-norm and will be denoted |H|/,2(f7)- 

We aim to prove error estimates compatible with the standard analysis of the finite element 
method; full details on the setting are available in a number of textbooks, e.g., [3l HJ 13 SO] For 
linear, Lagrange interpolants, the optimal error estimate that we seek has the form 

— < Cdiam(f]) |ii|^2(Q) , Vu G (2) 

where / is the interpolation operator lu := 'u(vj)Aj(x) with the summation taken over the 
element vertices. 

Since we consider only invariant (E(3]) generalized barycentric coordinates, estimate ^ only 
needs to be shown for domains of diameter one as passing simple scaling and rotation operations 
through the Sobolev norms yields the factor of diam(r2) for elements of any size. More formally, 
assuming the estimate ([2]) holds for all diameter one domains, the scaling argument follows in a 
similar fashion as seen in Remark 12.51 Let be a diameter one domain and ut the scaled function 
defined on TQ where T is a uniform scaling to a different diameter. The estimate is established 
by changing integration variables to a uniform domain (where ut{T{x)) = u(x)), applying the 
diameter one result, and scaling back: 

\UT - IutIhi^tq) = diam(™)'^-2 \u - 

< C diam(TO)'^~^ 1^1/^2(0) ~ ^ diam(Tr2)^ I^T|/f2('pn) • 

The final equality has an additional power of diam(rr2)^ compared to the equation from Remark l2.5l 
since it involves the iJ^-norm and the chain rule applies. The above inequality only addresses the 
iif^-seminorm, but the remaining lower order component of the H^-norm (the L^-norm) follows a 
very similar argument yielding a larger power of diam(Tr2). 

Using barycentric coordinates satisfying B3 under geometric restrictions GUJand G|5j it is suffi- 
cient to bound the H^-novm of the barycentric coordinates. 
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Figure 3: Notation for Corollary 13. 1[ 



Lemma 2.7 ([E]). For convex, diameter one domains satisfying CJI] o,nd CJ3 the optimal error 
estimate (0j holds whenever there exists a constant C\ such that 

||Aj||^i(f^) < Cx- (3) 

Lemma 12.71 is essentially the standard application of the Bramble-Hilbert lemma [2] in the 
analysis of the finite element method. While simplicial meshes only require a single estimate over 
the reference element, generalized barycentric coordinates need uniform estimates over all convex 
domains. Fortunately, the Bramble-Hilbert estimates can be shown uniformly over the set of unit 
diameter convex sets [351 E] a-nd thus the standard techniques apply. For a complete discussion of 
the framework and details of Lemma [2. 71 we refer the reader to [16]. Recalling that Gl5] follows from 
Gd] and G|2l the remainder of the paper is dedicated to verifying ([3]) for the mean value coordinates 
under G(l] and G|2] for the class of domains with diameter one. 



3 Preliminary Estimates 

First, we consider a simple fact about the constant /i* in Proposition 12.21 the sufficiently small size 
such that any ball of radius /i* does not intersect two non- adjacent edges of J7. 

Corollary 3.1. Under and CJH /i* < |vj — Vj„i|/2 for all i. 

Proof. Suppose the bound fails for some i. Then the ball -B ((vj + Vi_i)/2, /i*) intersects three 
edges of the polygon contradicting Proposition 12. 2| see Figure [3l □ 

The next two results apply Proposition 12.21 to show that rj(x) is small for at most one index i 
and ai(x) is large (i.e, near vr) for at most one index i. 

Corollary 3.2. Under and CJH if ri{-x.) < /i* then rj(x) > /i* for all j / i. 

Proof. Suppose ri{x) < h^. Then B(x,h*) intersects the two segments which meet at Vj. If 
rj(x) < then B(x, h^) would also intersect the two segments which meet at Vj and thus S(x, /i*) 
would intersect a total of at least three segments contradicting Proposition 12.21 □ 

In Proposition 13.31 we show that under our geometric restrictions, at most one angle ai(x) can 
be large for a given x. 
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Figure 4: Notation for Proposition \'S.'S\ Case 2. 



Proposition 3.3. Under (Jl\and (J^ ifai{x) > a* := max ^vr — /5=f/2, 2 arctan ^7^^^ then aj{x) < 
a* for all j ^ i- 

Proof. Fix X and suppose that aj(x) > a*. 

Case 1 : j £ {i — l,i + 1}. For the j = i — 1 case, consider the quad with vertices x, Vj_i, Vj, and 
Vj_|_i (see Figure m right). By condition G|3]and the fact that the sum of angles in a quad is 2tt, we 
have 

ai(x) + aj(x) + /?* < aj(x) + aj(x) + fii < 27t. 

Rewriting, we have that aj(x) < 27r — — aj(x). This estimate for aj(x) also holds when 
J = i + 1 by a similar argument for the quad with vertices x, Vj, Vj+i, and Vj+2- By hypothesis, 
Q;j(x) > vr — /3*/2, whence aj(x) < vr — /3*/2 < a*. 

Case 2 : j ^ {i— 1, i + 1}. Divide the triangle Av,Vj+ix into two right triangles as shown in Figured 
For the right triangle containing the vertex furthest from x, we adopt the notation of the figure: 
let Oi be the angle incident to x and let hi and ii to be the lengths of the two sides depicted. By 
choosing the furthest vertex, 9i > aj(x)/2. Since tan^j = ii/hi, 

hi = — - — ^ ^ /i* . 

tan 6i tan (oj (x) /2) 

The final inequality above results from our assumption that aj(x) > 2 arctan So i?(x, h^,) must 
intersect the segment between Vj and Vj+i. By Proposition 12. 2| B{x,h^:) cannot intersect the 
segment between Vj and Vj+i (because that segment is not incident to Vj or Vj+i). 

Now define 6j, £j and hj in a similar fashion to 9i, ii, and hi, except corresponding to the 
segment between Vj and Vj+i. Since B(x, h^) doesn't intersect vyv^IjlT, hj > /i*. Then aj(x) < 29 j 
and 

tanOj = -L <—. 

j * 

Thus aj(x) < 26lj < 2 arctan (^-^^ < a*. □ 

The next two results prove some intuitive notions about the size of aj(x) when x is near the 
boundary of The first (Proposition 13. 4( ) says that a 'big' Oi value and 'small' rj value can only 
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Figure 5: Notation for Proposition 13.51 

occur simultaneously if Vj and Vj are identical or adjacent. The second (Proposition 13. 5p shows 
that if X is close to a vertex, the two aj angles defined by the vertex have a 'large' sum. 

Proposition 3.4. Under and if Oii{x) > a* and rj(x) < /i* then j G + 1}. 

Proof. As we saw in Proposition [3i3l if ai(x) > a* then -B(x, h^,) intersects the line segment between 
Vj and Vj+i. Thus Proposition 12.21 ensures that B{x,h^:) cannot contain Vj for j ^ + 1}. □ 

Proposition 3.5. Under (J^and C@ if ri{x.) < /i* then ai_i(x) +Qj(x) > 27r/3. 

Proof. Define := Zxvj_iVj, Q := Zxvj+iVj, and recall /3j := Zvj_iVjVj+i; see Figure [5l By 
Corollary 13.11 we have ri < h* < |vj_i — Vi| /2. By the law of sines, 



sin^j 



Vi_l 



■ sinai_i < 



sm ai-i 



Hence, < arcsin y — — Similarly C,i < ai/2. Summing the interior angles of the 

quadrilateral with vertices x, v^^i, Vj, and Vj+i gives 



ai-i + ai + ^i + Ci + /3i 



2tt. 



Since /3j < vr, we have a^-i + + + Ci ^ Applying the inequalities on and Q gives the 
result. □ 

Thus far, we have given bounds on the size of angles for a fixed x value. In the next section, 
we will also need estimates of how fast Q!i(x) is changing, i.e., estimates of |Vaj(x)|. The next 
proposition provides an estimate on this term. 



Proposition 3.6. |Vai(x)| < 



1 



+ 



1 



ri(x) ri+i(x)' 



Proof. Without loss of generality, let Vj = (0,0) and let Vj+i = (d, 0). We will establish this 
estimate for any d. Also let {x,y) := x. Define Oi, rji as shown in Figure [6] so that ai(x) = 9i + i]i 



with tan = | and tan rji 



Differentiating 9i with respect to x and y, we find that 



1 + {x/yY 





1 


y 


_ -x/y"^ _ 


_j_ y2 


— X 
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Since ri{x.f = + y^, it follows that |V6'i| = + = Similarly, \V7]i\ = 

We note that these estimates on 6i and rji are independent of the edge length d: they only depend 
on the locations of Vj and Vj+i, respectively. As Vaj(x) = V9i + Vr/j, the triangle inequality 
completes the proof. □ 

Since rj increases radially from Vj, we also have a simple bound on the gradient of r^. 
Proposition 3.7. Vrj(x) = j -j- and hence |Vrj(x)| = 1. 

Our final result of this section is a conservative uniform lower bound on the sum of the weights 
Wi at an arbitrary point x. This ensures that the denominator of the mean value coordinates {A.j} 
never approaches zero. 

Proposition 3.8. X]"=iW'i(x) > 27r. 

Proof. Since our domain has diameter 1 (see Remark 12. 5p . we have rj(x) < 1. Thus 

ran ran 

^Wi{x) > 2 J];tan(a,(x)/2) > 2^Qi(x)/2 > ^ai(x) = Itt. 

i=l i=l i=l 1=1 

□ 



4 Main Theorem 

Our main result, Theorem 14. 3( is a uniform bound on the norm of the gradient of the mean value 
coordinate functions under Gl and G2. The proof works by writing 



2 



where A'^i and N2 are given in terms of {tj} and {rj}. The summands in Ni and A^2 are bounded 
by constant multiples of (X^i^j)^! as shown in Lemma l4.ll and Lemma 14.21 respectively. 
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Case 2 

Figure 7: Division into cases for Lemma |4.1[ 



Lemma 4.1. Under conditions and and for a ^ b, there is a constant Ci such that 

2 



(ta-l(x) + ta(x))(t6„i(x) + tfe(x)) 



Vr„(x) 



ra(x)2r6(x) 



(4) 



/or aZZ X G . 

Proof. Fix X E ri. The argument is separated into two cases based on the distance from x to v^; 
see Figure [71 We will make use of the facts that for any index i, 

''i(x) < 1 and |Vrj(x)| = 1, 

by the diameter 1 domain assumption and Proposition 13.71 respectively. For readability, we omit 
the dependencies on x from the explanations. 

Case 1 . ra(x) > /i*, i.e. x away from v^. 
Since 



{ta-l + ta){th^i + tfo) = {Wara){Whrb) < [Wa + Wb) rarb, 



it follows that 



[ta~l + ta)[tb-l + tb)^ — 

rirb 



< {Wa + Wby 



Wrn 



< 



Case 2 . ra(x) < /i*, i.e. x close to Vq. 
By Corollary 13.21 and Propositions 13.31 13.41 and 13.51 we conclude: 

rfe(x) > /i* 

min(aa_i,aa) < a*, 

Oj < a*, for i ^ {a — l,a}, and 

max(aa_i(x),aa(x)) > tt/3. 



(5) 
(6) 
(7) 
(8) 
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Let m := max (0^-1 (x), aa(x)). We break into subcases based on the size of m relative to a* (as 
defined in Proposition I3.3p . 
Subcase 2a : m > a* . 

By ([6]), ([7|), and the subcase hypothesis, we have Ui < m for any i. Since m = a^-i or m = Oa, 
we have tan(m/2) < ta^i + ta- Hence 

h-i +h<2 tan(m/2) < 2 (ta-i + ta) • 

Using this and ([5|), we conclude that 



< 2 



^ 2 

/IT"- 



Subcase 2b: m < a* . 



By ([6|), d?]), and the subcase hypothesis, we have Oi < a* for any i. By dSj), tan(7r/6) < tan(m/2) 
and hence tan(7r/6) < ta^i + ta- Putting these facts together, we have that 



tb-i + tb< 2tan(a72) < 2tan(a72) 



tan(7r/6) 



Using this and we conclude that 



{ta-l + ta){tb-l + tb)—^ — 



< 



(ta-l+ta)^ 2tan(QV2) 



tan(7r/6) 



2tan(a*/2) 2 



/i*tan(7r/6) " 



In both subcases, the observation that < ( ^ ■ Wj{x) j completes the result 



□ 



Lemma 4.2. Under conditions CJIl o'^c^ o'^c^ /o?" * 7^ J 0?^-^ cl ¥^ b, there is a constant C2 such 
that 



Vti{^)tj{^) 



ra(x)rfe(x) 



(9) 



/or all X & fl. 



Proof. Fix X £ ri. For readability, we will often omit the dependencies on x from the explanations. 
By Corollary 13.21 at least one of {ra,rb} is bigger than h^. Without loss of generality, assume that 
i^a ^ fb- Similarly, by Proposition 13.31 a-t least one of {ai,aj} is smaller than a*. 

Since the left side of ([9]) is not symmetric in i and j, we must break into a number of cases 
based on both the comparisons of Oi and Oj to a* and the comparison of to /i*. The regions 
where each case holds are shown in Figs. [8]and[9l 

Note that in many of the cases, we will make use of the fact that a* > 7r/2. This is confirmed 
by starting with the trivial bound /?* < vr and then deriving 7r/2 < vr — /3^,/2 < a*. We will also 
frequently make use of the following bound on Vti(x). Observe that 



\Vti{x) 



V tan 



Q!i(x) 



^|VQi(x)| sec^ 



|VQ;i(x) 



2cos2 (ai(x)/2)' 



Using the bound on |Vaj| from Proposition 13.61 we get the bound 



\Vti\ < 



1 

— + 



2cos2(ai/2) \ri rj+i/ ' 



(10) 
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3 and 4 ,< 3 and 4 




Figure 8: The proof of Lemma [4. 2 1 is broken into numbered cases according to where x hes relative 
to vertex Vq. The middle arc is the radius h^, ball around v^. Inside the other arcs either aa-\ or 
Q.a is bigger than a*. 




Figure 9: Example configurations of vertices for different subcases in Lemma 14.21 



Finally, we will make use of an additional index k defined by 

k := argmin{ri,ri+i}, 

i.e. Tfc is the shorter of rj and rj_|_i. 
Case 1 . aj(x) > a* and ra(x) < /i*. 
We immediately have aj < a* and > h^, and hence 



(11) 



< 



tan(QV2) 



(12) 



Since tt > > a* > tt/2, we have sin(aj/2) > sin(7r/4) = l/V^ or, equivalently, 1 < 2sin^(aj/2) 
This fact, along with (fTO|) . gives us the bound 



iVtil < 



1 _^ 1 \ sin2(ai/2) ^ 2 \ „ 



2cos2(ai/2) ^rj ~^ r^+i J ^ cos'^{ai/2) ^ra) r„ 



By Proposition 13.41 i S {a — 1, a} meaning ti < ta + ta-i = WaTa- Thus |Vtj| < 2w'^ra. Combining 
this estimate with (|12p . we have 



ran 



2tan(a*/2) n 
< r — —wi. 
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Case 2 . aj{x) > a* and ra(x) < /i*. 

Proposition 13.31 and Corollary 13.21 imply that < a* and ri, > h^, > ra- Since < < a* < tt, 
we have 1 > cos(Qi/2) > cos(a*/2) > 0. Combining these facts with (jlOp gives 



iVtil < 



1 



1 1 

+ 



< 



1 



2cos^(aj/2) \rj ri^i J 2cos2(a*/2) raCos2(a*/2) 



Since a* > 7r/2, we have tj > tan(a*/2) > 1. By Proposition 13.41 j £ {a — l,a}, allowing the 
bound tj < t'j < {ta~i + ta)^ = w'^r'^- Putting all this together, we have that 



Vtitj 



ran 



< 



< 



wi. 



raCOS^(a*/2) r^h* cos^ {a* / 2)h^ " 



Case 3 . aj(x) > a* and ra(x) > h^. 

As in Case 1, we have 1 < 2sin^(aj/2) so that 



Vtitj 



Tan 



< 



tan(a*/2) 



2cos2(ai/2) \ri r. 



1 1 

+ 



tan(aV2) 2 ^ 



Since < diam(r2) = 1, |rfc|^ < 1. As k £ + 1}, we have tf/\rk\ < {tk-i + tk) /Ifkl"^ < wl- 



Thus, 



2tan(a72) o 



/i2 



Case 4 . aj(x) > a* and ra(x) > /i^,. 

By the same arguments as in Case 2, we have 



ran 



< 



2cos2(a*/2)/i2 



ri ri+i/ cos2(a*/2)/i2 



Subcase 4a . > /i*. Since a, > a*, t,- < < lo^r^ < i«f . Thus 



Vtitj 



rart 



< 



1 



cos2(a*/2)/i3 



It; 



Subcase 4b . < /i*. Proposition 13.41 implies k G {j,j + 1} and hence tj < t^ < (tfc-i + ifc)^ 



wlrl- Thus, 



< 



jr, < 



1 



cos2(a*/2)/i2 ^ ^ - cos2(a*/2)/i2 ^' 
Case 5 . ai < a* , aj < a*, and Va < h^^. 

As before, we begin recalling that Corollary 13.21 implies < rk- By Proposition 13.51 ta-i +ta > 
2tan(7r/6) > 1. (Note: tan is a convex function function on (0,7r/2) and thus the smallest value 
occurs when Ua-i ~ ^ 7r/6.) Then using (fTO]l . we estimate 



Vtitj 



rarb 



< 



t&n{a*/2) /I 1 \ tan(a*/2) 1 ^ tan(aV2) ^ 
2 cos2(a*/2)/i*r(j \rj J'j+i / ~ cos2(a*/2)/i^ ~ cos2(a*/2)/i^, " 
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Case 6 . ai < a* , aj < a*, and ra> h^,. 

First, following similar estimates to previous cases yields, 



Vtitj 



ran 



< 



tan(a*/2) /I 1 \ tan(aV2) 1 
2cos2(a*/2)/i^ \ri rj+i/ ~ cos2(a*/2)/i2 



Subcase 6a . > h^. By Proposition 13.81 we have that 



Vtitj 



ran 



^ 2tan(a*/2) ^ 2tan(a*/2) 



cos2(aV2)/i3 - 47r2cos2(a*/2)/i3 



Subcase 6b . < /i*. By Proposition 13.51 t^-i + ifc > 2tan(7r/6) > 1. Thus, 



Vtitj 



rarb 



tan(a*/2) 1 tan(aV2) tan(a*/2) 2 

- cos2(aV2)/i2 ■ ^ - cos-^ia* /2)hf^ " cos2(aV2)/i2 "^^^^ 



In each case/subcase, the desired estimate holds. Taking the maximum constant over each case 
completes the proof. □ 

Theorem 4.3. Under conditions (J^ and tJH there exists a constant C such that 

|VAi| < C 

Proof. For readability, we omit the dependencies on x from the explanations. By the quotient rule, 
the gradient of a weight function Wk can be expressed as 



Vtfc-i + Vtfc tk-i + tk 



Similarly, the gradient of Aj can be expressed as 



VA,; 



(13) 



(14) 



Plugging ()13p into ()14p , we partition the summands of the numerator according to whether or not 

2 

I -^1 + -^2 1 where 



they involve some Vr^ factor. We thus write (^^jWj ] \VXi 



Ni = Y{ti-i+ti){tj_i + tj 



Vrj Vri 



1 2 
r^ri rfri 



and 



n ^ 

N2 = Y1 ii^ti-i + Vti){tj.i + tj) - {ti_i + ti){Vtj.i + Vtj)] . 

To bound |A'^i|, note that the i = j terms cancel and there are at most 2n* terms in the summation. 
Thus Lemma |4. II applies and we have 



|iVi| <2n*Ci [Y^wA . 
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Figure 10: A simple computational example is given for a mesh of "degenerate octagons", i.e., 
squares with mid-side nodes (left). Basis functions corresponding to a corner node (center) and a 
mid-side node (right) are shown. 

To bound |A'2|, note that it can be expanded into at most 8n* terms of the form 

ran ' 

The terms with k = I or a = b cancel each other out meaning Lemma 14.21 applies. Thus, 



where here \Q\ denote the area of ft. Thus, by Lemma [2. 71 Theorem 14 . 3 1 guarantees that the optimal 
interpolation error estimate ([2|) holds. 

5 Numerical Example and Concluding Remarks 

By bounding gradients of the mean value coordinates uniformly over the class of polygons, we have 
formally justified one of the key motivations for the use of the coordinates. Moreover, this bound 
is the essential ingredient in the optimal interpolation error estimate. We briefly demonstrate that 
our interpolation result translates to standard convergence of a finite element method using a mean 
value interplant operator. To demonstrate success of the mean value basis in the presence of large 
interior angles, a mesh is constructed of "degenerate octagons", squares with additional nodes in 
the middle of each side; see Figure [TOl With a basis of mean value coordinates, we solve Poisson's 
equation with Dirichlet boundary conditions corresponding to the solution u{x,y) = sin(x)e^. As 
shown in Figure [TTl the expected convergence rate from our theoretical analysis ([2]) is observed. 




Putting these together, we have 



VAi(x)| < 2Cin* + 8C2n* 



which is the desired bound. 



□ 



Finally, note that Theorem 14.31 implies ([3]): for a diameter one domain. 
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\\U- Uh 


L2 






n 


error 


rate 


error 


rate 


2 


3.35e-3 




7.56e-2 




4 


8.67e-4 


2.03 


3.60e-2 


1.07 


8 


2.18e-4 


1.99 


1.76e-2 


1.03 


16 


5.50e-5 


1.99 


8.73e-3 


1.01 


32 


1.38C-5 


1.99 


4.35e-3 


1.00 


64 


3.47e-6 


1.99 


2.17e-3 


1.00 


128 


8.69e-7 


2.00 


1.09e-3 


1.00 



Figure 11: Uniform refinement of a sequence of degenerate octagonal meshes yields the expected 
convergence rate using mean value basis functions. Meshes of elements are shown for 7i = 2 
(left) and n = 4 (center). Tabulated results (right) for the solution of Poisson's equation with 
Dirichlet boundary conditions demonstrate second-order convergence in the L^-norm and first order 
convergence in the i/^-seminorm. 



namely, linear convergence in the H^-norm. The quadratic convergence in the L^-norm is also 
expected from the Aubin-Nitsche lemma; see e.g. (3j. 

Another advantage of mean value coordinates is the fact that the formula can be evaluated for 
non-convex polygons, while some other coordinates (e.g., Wachspress) are not defined. While mean 
value coordinates can become negative for certain non-convex polygons (especially in the presence 
of interior angles near 27r), the interpolants are satisfactory in some applications [18]. To get the 
gradient bound in Theorem 14. 3( convexity is only used in a few places. Specifically, Proposition 12.21 
is not true for general non-convex sets. Instead, analysis in this setting should be restricted to 
the class of non-convex polygons for which a constant /i* > exists such that i?(x, /i*) does not 
intersect three polygon edges. Additionally, Proposition 13.31 fails: a point may form large angles 
with two adjacent edges when the edges form a large (near 27r) interior angle. While pinning 
down precise geometric restrictions for bounded gradients on non-convex polygons becomes overly 
complex, our analysis does give some intuition as to why mean value coordinates succeed in many 
common applications involving non-convex regions. 

Finally, mean value coordinates can be defined for 3D simplicial polytopes [15] (in addition 
to a Wachspress- like construction [21]). While we expect that a similar analysis of interpolation 
properties can be performed in this setting, there are two primary obstacles. First, precise 3D 
geometric restrictions must be posed which can become rather complex; dihedral angles must be 
considered in addition to the quality of all simplicial facets. Naive hypotheses can lead to an overly 
restrictive setting. Second, the 3D analysis will involve many more cases than the already involved 
2D analysis. A better approach may be to identify new generalizations that simplify the existing 
proof before extending the results to 3D. 
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